SPECTRAL ANALYSIS OF A CLASS OF HERMITIAN 
JACOBI MATRICES IN A CRITICAL (DOUBLE ROOT) 

HYPERBOLIC CASE. 



SERGUEI NABOKO AND SERGEY SIMONOV 

Abstract. We consider a class of Jacobi matrices with periodi- 
cally modulated diagonal in a critical hyperbolic ("double root") 
situation. For the model with "non-smooth" matrix entries we ob- 
tain the asymptotics of generalized eigenvectors and analyze the 
spectrum. In addition, we reformulate a very helpful theorem from 
[9] in its full generality in order to serve the needs of our method. 



1. Introduction 

In the present paper, we consider a class of Jacobi matrices with peri- 
odically modulated growing diagonal giving an example of the " double 
root" problem. This means a critical situation for asymptotics of the 
generalized eigenvectors related to the matrix. Such a situation partic- 
ularly arises in the spectral phase transition phenomenon. If the matrix 
depends on some parameters, the decomposition of its spectrum into 
different types (absolutely continuous, singular continuous, pure point, 
discrete) may be independent of these parameters. But if the structure 
of this decomposition changes under a variation of the parameters, a 
spectral phase transition occurs. If this change happens by a jump, 
such a phenomenon is called a spectral phase transition of the first 
type, whereas if this change is smooth with the change of the parame- 
ters when they move across some hyper-surface in a space, it is called 
a spectral phase transition of the second type. The transition in types 
of spectrum is tightly related to the change of form of asymptotics of 
generalized eigenvectors due to subordinacy theory (Gilbert-Pearson 
[6]), which was generalized to the case of Jacobi matrices in [12] . 
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Last decade the spectral analysis of Jacobi matrices attracted the 
attention of many specialists in operator theory and mathematical 
physics. 

We consider the Jacobi matrix J with diagonal entries 

for odd values of n 
for even values of n 

(where a and h are real parameters a G (|; 1) and 6 7^ 0) and off- 
diagonal entries (weights) 

a„ = ra". 

The model demonstrates the situation of the spectral phase transition 
of the first order corresponding to the "moment of transition" exactly. 
We are interested in asymptotics of generalized eigenvectors, i.e., solu- 
tions of the spectral recurrence relation 

(1) a„_iM„_i + hnUn + a„Mn+i = n>2. 

This model has been studied in |5] and in particular the following result 
was obtained: for A < there are two solutions and u~ of ^ with 
the following asymptotics as n ^ 00: 



4n ~ (-l)"n-f exp (±^^^ 

±\/5 (1 - f ) (-l)"^-"^ exp {±^^^ 



The problem of determining the asymptotics for A > was stated, 
and in the present paper, we show that the answer has the same form. 
However this interesting question is not the main concern of this paper. 
The principal difficulty in our analysis is that the situation is "critical 
hyperbolic" unlike the "critical elliptic" situation in [5], which roughly 
means that exponents in the answer grow and decay if A > and 
oscillate if A < 0. 

Let us explain the problem of the critical situation in more detail. 
Whenever one deals with the three-term recurrence relation ([T]), it is 
often useful to write it in the vector form introducing the sequence 

'^n '■= ( ^"^^ ) and the transfer-matrix i?„ := ( _a„-i \-b„ )• So 

V " / \ an an / 

dl]) is equivalent to the discrete linear system in C^: 

(2) l^n+l = Bnltn, U > 2. 

The solution for such a system is obtained by taking the chronological 

product of transfer matrices, (11-^^) "^2 = -Bn ■ ••• ■ -B2 ^^2- The 

\fe=2 / 

analysis becomes much easier if the matrix of the system is in some 
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sense smooth-in-n (say, has a limit and asymptotic expansion in inverse 
powers of n as n — 7- oo). In our case the transfer-matrix is not smooth 
in this sense because the coefficients of the spectral equation "jump" 
all the time. This is the first (simple) problem that we face and it may 
be solved by taking the product of two consecutive transfer-matrices, 
introducing the new linear system with the sequence of the coefficient 
matrices 

(3) Mn := B2nB2n-l 

which turns out to be smooth-in-n. It is easy to see that the sequence 
M„ has a limit M := lim M„ with detM = 1. The are three possibil- 

n—^oo 

ities for the eigenvalues of M: they can be 

• unimodular and complex conjugate (the elliptic situation), 

• real and different (the hyperbolic situation; hence one of them 
is greater than 1 and another less than 1 in absolute value), or 

• coincide and equal 1 or —1 (the critical situation=the double 
root case). 

In the hyperbohc situation solutions are supposed to grow or decay, 
in the elliptic situation they are supposed to oscillate having similar 
behavior of their norms The numerous variants of analogues 

to Levinson Theorem (also known as Bernzaid-Lutz Theorem [4j) for 
differential linear systems [3] can be applied in these two cases, cf. [9], 

m- 

In the critical situation one can also distinguish the " critical elliptic" 
and the "critical hyperbolic" cases. This separation depends on the 
lower orders behavior of M„ as n — )■ oo, namely on the asymptotic sign 
of the discriminant (discrM^ := (trM„)^ — 4detM„) of matrices M„. 
In the critical case, the matrix M is similar to the Jordan block, with 
powers 

1 J V 1 
and the difficulty here is not that the powers of the matrix grow with 
n (as they do in the hyperbolic situation), but the fact that the large 
entry is off-diagonal. This mixes upper and lower components of the 
solution, which makes the system unstable and sensitive to small per- 
turbations. The problem of the "critical hyperbolic" situation was 
considered for smooth matrix elements in [10]. In the present paper, 
considering a model of Jacobi matrix with "oscillating" diagonal, we 
intend to simplify greatly the scheme of successive transformations of 
the matrix system which was used there by making it more general 
and transparent (Section [3]). This approach differs from that of ^ for 
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"critical elliptic" situation, because in our case we have to deal with 
growing exponents (see Step 3 on page[TD]). Moreover, we consider in 
the appendix a theorem from [9] which is necessary for the final step 
of our method. The problem is that the original formulation of this 
theorem states asymptotics of only one (principal) solution. Although 
it is surely the most difficult part of the problem considered, we want 
to state the result in its full generality. 

The method of the present paper works for a much wider class of 
Jacobi matrices. However our goal is to present a simple formulation 
of the method and show by means of an example how it works. 

A similar problem (the critical hyperbolic situation) was considered 
in [7] for a " smooth" model where the author used a completely differ- 
ent method related to Kelley's paper pT| . 



2. Preliminaries 

As usual, the operator J in the Hilbert space /^(N) is first defined 
(as J) on the linear set of vectors which have only a finite number of 
non-zero components, //^(N), by the rule 

{Ju)i = OiMi + biU2 

{Ju)n = bn-lUn^l + ttnUn + &nMn+l, for U > 2. 

Then its closure J = is a self-adjoint operator provided the Carleman 

oo 

condition J2 ~ — +^ i^J is satisfied. In the standard basis {en}'^=i 

n=0 

(where e„ is the vector with all the components zeros except the n-th) 
the operator J admits the following matrix representation: 



/ h 


ai 


■ 


■ \ 


ai 











a2 















We start with the system 



( ]=mJ 

\ M2n+1 J \ U2n-l J 
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see 



(j3])). One can obtain a solution of this system directly, by 



taking a product of matrices M„: ( 
trices M„ are: 



n Mk 

k=2 



U2 

Us 



Ma- 



B2nB2n~l 



-1 



+ 



A 





-1 



-1 y (2n)" V -1 

having smooth-in-n asymptotics as n — i- oo. 
The eigenvalues of the limit matrix M := 



+ 7^1 + { — 



2n 

-1 -b 
-1 



n 



2a 



coincide, so 



the situation is critical (the double root case). By an easy calculation, 
one can see that discriminants of matrices M„ equal discrM„ 



4fcA 

i2ny 



+ 



O (^), hence indeed A < corresponds to the elhptic situation and 
A > corresponds to the hyperbolic situation (both critical). 

Our method is based upon a sequence of transformations which are 
determined by some anzats. Similar consideration for the anzats con- 
structing may be also found in [TO] . 



Remark 1. In what follows, we use transformations that are in fact 
discrete analogues to variation of parameters method transformations. 
Whenever one deals with the product of matrices, say An, one can as 
well consider matrices Cn = T~_^^AnTn (where the matrix sequence Tn 
is chosen somehow; we will call matrices T„ affinity-like) . Due to the 
cancellation of intermediate terms the product of matrices equals 



"2 



"2 



•n,2 



n=ni n=ni \n=ni / 

So the study of the linear difference system with coefficient matrices An 
can be completely reduced to the study of the linear difference system 
with coefficient matrices Cn- 

3. Calculation of the asymptotics in the hyperbolic case. 

In this section, we proceed through several transformations in order 
to simplify the problem and finally obtain the system which can be 
treated with Janas-Moszynski Theorem (which is properly adjusted in 
the appendix). So we divide this section into four steps. In fact we 
have already made a "zero" step in the first section, which is reduction 
to the smooth matrix system by taking the product of two transfer- 
matrices. But we do not include this step into Section |3] because the 
double root problem actually arises only at this stage. Moreover after 
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grouping the transfer-matrices by pairs one does not need to perform 
any inverse transformation in order to obtain the answer. 

Step 1: Reduction of the meaningful part of M„ to the transfer- 
matrix form. 

We write the spectral equation in matrix form to settle the problem 
of periodically modulated coefficients by grouping transfer-matrices in 
pairs. But now it is simpler to consider a smooth three-term recurrence 
relation which is equivalent to the system. The transfer-matrix corre- 
sponding to the three-term recurrence relation should have entries 
and 1 in the upper row. In this step, we find the proper transformation 
which makes the coefficient matrix of the system resemble a transfer- 
matrix. This transformation is generated by the matrix sequence of 
the form 



Tn := (-ir [( I 
The transformation 



A 



2b 





















(5) 



T 



n+1 



MnT- 



gives 



Nn 



1 
-1 2 



b\ 
(2n)" 




1 



+ 



a 



n 



+ 



n 



The meaningful part (or the part which we expect to be meaningful - 
the sum of first few terms) is now a matrix with the upper row of the 
form 0, 1. The exact form of matrices T„ can be determined from this 
requirement (but the answer is not unique of course). Indeed, try to 
find Tn in the form 



(6) 



T 

■J- T). 



1 -b 

1 



^ yd) I ^2-(2) 
(2n)" 2n 



with unknown matrices T^^^ and T'-^^ independent of n. Denote T : = 
J ^ j which is chosen in order to satisfy 



-T 







T 



-1 



1 
-1 2 



=: N. 



N T 



By substitution of T„ in the form ^ into the relation T„+iM„ 
looking at the terms of orders 4r and - one obtains the following linear 
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By stars we denote matrix entries which are allowed to be non-zero. 
Therefore the problem may be reduced to the following one. To prove 
for the commutator equation 



that for any given values fi and /2 there exist unique values xi and X2 
and 2x2 matrix X (obviously not unique) which satisfy the equation. 

Multiplying the equality on the left by the matrix ^ ^ ^ 

on the right by ^ J ^ ^ (its inverse) we obtain another form of the 
commutator equation: 



-y-i yi-yi \ ^ ( /i + /2 /2 

ys J \ Xi + X2 ~ fl - f2 X2 - f2 

We denoted F := I \ _ / 1 ^ x f ! ° 1 . It follows 



1/3 ^4 / ■ \ -I I J V ^ ^ 

now that X2 = —fi and xi = 2/i + f2 are uniquely determined by fi 



and /2, and the matrix Y exists and is unique up to ciJ + C2 



1 




with any Ci and C2- 

Therefore one can reduce the original system (jlj) to the new one, 
where M„ are replaced by iV„. Note that the exact form of the matrices 
T^^^ and T^'^^ is not essential for the result of transformation. We need 
only to prove the existence of such matrices. 

Step 2: Reduction of the main term to the identity matrix. 

The structure of the main part of the system allows to write 



^^"=1 -F°(«) -Fdn) ^^''^ ' 



n2" 
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where the error term is a 2 x 2 matrix which norm is 0(?2~^") and the 

bX \ 



matrix entries Fi{n) and Fi{n) are (denote B := \^ 

Fi n =-2- — + -, F,{n) = l--. 

n n 

Our goal is to use the fact that the system 

Wn+i ) \ -F2{n) -Fi{n) J \ Wn 

where the remainder is omitted, is equivalent to the three-term recur- 
rence relation 

(7) Un+l + Fi{n)un + F2{n)Un-l = 0. 

The latter has two "approximate solutions" of the form 



with 7 = — ^,5 = 1 — ^,^4=^. This means that 

(8) + F^in)z^ + F,{n)zt, = 0{n-^^)z^ 

and can be verified by a direct calculation. The form of these "ap- 
proximate solutions" can be obtained, for instance, by analogy with 
the WKB method (see P2] for that type of argument). Having this 
structure of the "solution" one can determine unknown values of 7, 5 
and A from the condition of cancellation of all decreasing terms in (jHI) 
up to the order 0(?2~^"). Equation ([8]) implies that 



(9) 



\ // 1 \,/0 



4+1 ; " VV -^2(n) -F,{n) ; + V 0{n-^-) ) ) \ z. 

= (iV„ + 0(n-^")) 



± 



"71-1 



zt 



It is useful now to write the last vector equality in a matrix form. 
Denote 



^n-l ^n-1 



5„ = 

then combining ^ for both signs one has 

5n+i = (iVn + 0(n-2-))5„ 

and hence 

This actually follows only from the fact that "approximate so- 

lutions" of the recurrence relation ([7]), i.e., from (|8]). 
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Substituting the expression for = rCe^'^'^ into the term 

and using the fact that det Sn+i ~ 2y45n~" as n — )• oo (for calculations, 
cf. [5]; one can also find this rate of decay from the fact of modified 
wronskian persistence), one obtains 

(10) SM^I.'fOW 



^n+l^^n^n ^ ^3^/2 ^ e-2^"'0(l) 0(1) 

These calculations enable us to reduce the original system to the 
system with the coefficient matrix Kn, for which the "main term" is 
identity matrix. The problem is that it contains exponentially increas- 
ing anti-diagonal terms. In the next step we show how to overcome 
this difficulty due to the symmetric cancellation of both anti-diagonal 
entries. 

Remark 2. The above calculations demonstrate a simple "geometrical" 
approach based on reduction to the tree-term recurrence relation ([7]). 
Another approach (a simplification of calculations from [10\) may be 
rather straightforward: the "geometrical construction" is replaced by an 
explicit calculation. Since det Sn+i = -2^-2^+1 — -2^-2^+1, the substitution 
of matrices Sn gives: 



^^ni^n+l + Pl^n + -^22:„_i) -Z^Z^^i - Z~{FiZ^ + -F2Z„_J 

After adding and subtracting terms z^z^^-^ in the upper-left entry and 
-2~-2^+i in the lower-right entry for extracting the determinant, the last 
expression becomes the following one: 



X 



~^ni^n+l + -^l^n + -^2^„_l) -^,7(^„+l + -^l^n + -^2^„,_l) 



It is remarkable that the expressions z^_^_i+ Fiz^ + F2Z^_i appear in each 
matrix entry of the second term. Now substituting ([8]) and taking into 
consideration that det Sn+i ~ 2A6n~°', one obtains the same expression 
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as m ([TOD : 



detSn+A z-'0{n-'^) ztz-0{n-^'^) 

1+ - i 0(1) e^^-'Oil) 



3a/2 \^ 0(1) 0(1) 



The problem now is the growing exponent in the upper-right entry 
of the matrix, which leaves no chance to regard the remainder small in 
any sense. This exponent can be compensated by the decaying one in 
the lower-left entry as we show below. 

Step 3: Elimination of exponentially increasing off-diagonal terms. 

In order to produce the elimination we perform yet another trans- 
formation with affinity-like matrix 



n 



g2A„* 

1 



which yields 

(11) L„ := X-^iir„X„ = ^ Q ^ j+0 



3a/2 



n 



1 2A<5 , {2A5f n. , 



1 ) \n^'^l^ 

As a result of all the transformations, the original system is reduced 
to a new one with coefficient matrices L„ for which Janas-Moszynski 
Theorem |9j (Theorem |2l Appendix) is applicable. Take p„ = 

Vn = V = ^ Q ^ and let -R„ = 0(72^^"/^) be the matrix remain- 
der which belongs to for a > 2/3. The theorem asserts that the 
system with coefficient matrices L„ has a basis of solutions of the form 

^-2An^ ("^1 + o(l)) and + o(l) (we use notations "e^i = 



^2 



Remark 3. Note that we need to apply Theorem\^ only to ignore the 
remainder term of order 0(n~'^"/^). Without it, the system generated 
by matrix ( [TT]) obviously has two exact solutions e~^"^" and ~^2- 
However this is not entirely obvious, since the eigenvalues of the limit 
matrix lim L„ coincide. 

n— ^oo 

Step 4: Returning to the original system. 
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Returning to the system with coefficient matrices M„ we recall steps 
1, 2 and 3. We should take into account the affinity-like matrices that 
stay in front of the product of matrices L„. Since 

one has (for every A we choose a number no such that detT„ 7^ for 
n > Uq): 

U2n 
U2n+1 




\k=no 



no+l 



-1) 



n+1 



4(l + o(l)) z+(l + o(l)) 



X 

.k=no 

As one has from the previous step, the vector 

n Lk] X^^^Sl^^Tno i J is a linear combination of e~^'^'^''("^i + 

k=no J V ^"0+1 / 

0(1)) and ("^2 + 0(1)). Therefore the vector ( ) is a linear com- 

\ U2n+1 J 

^ ^ { 2±(l + o(l)) \ 
bination of ( — 1)" , 4a^± /.nx • So the following result holds 

V ±6i^(i + o(i)) ; 

true. 

Theorem 1. For any A > spectral equation ([1]) has a basis of solu- 
tions and u~ with asymptotics as n ^ 00 of even components 

utn ~ (-1) n 4 exp ± 




2" 1 - f 



and asymptotics of odd components 



A «\ , 3a / bX n 



4+1 ~ ±\/ £r ( 1 - ^ ) (-l)"n-^ exp ( ± 



2^6 V 2 / ' ' ' \ V 2" 1 - f 




Regarding the operator J this means (due to the subordinacy theory 
of Gilbert and Pearson [6] [I2]), that the spectrum of the operator on 
the positive semiaxis is of pure point type. Moreover, if A G (Jp{J), 
then the solution u~ is an eigenvector of J. 
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Note that the asymptotics of as n — > oo for A > formally 
coincide with the ones obtained in [5] for A < as it was mentioned 
above. 

4. Appendix: Janas-Moszynski Theorem revisited. 

Let us turn our attention to a theorem from [9], which was not for- 
mulated in its full generality there. By using a complicated technique, 
the existence of the principal ("smaller") solution for some system was 
proved. But the existence of the second ("larger") solution was not 
stated. In what follows we prove this fact. 

Remark 4. For the sequences {an}, in order to avoid any non-essential 
problems related to vanishing of elements of the sequence, we assume 

N 

the notation for the product Y[ o.n meaning Yl ^n- We also use the 

11=1 7i;a,i7^0 

notation and for matrix sequences, i.e., the sequence of matrices 
{Mn}'^=i belongs to 



• l\ iff E < oo, 

n=l 

oo 

. D\ iffY. ||M„+i-M„ 

n=l 

Consider the linear difference system in 



(12) ]={I + PnVn + R„ 

\ '-'n+l 

Suppose it is non-degenerate, i.e., det(/ + PnVn + Rn) 7^ for every n 
(this means that the system has two linearly independent solutions). 

Theorem 2. (see [H] where the the essential part of the theorem was 
proved). 
Let 

oo 

• Pn be a positive sequence such that J2 Pn = +oo, 

n=l 

• {Rn} ^ l^ be a 2 X 2 matrix sequence, 

• {Vn} E be a real 2x2 matrix sequence with discr ( lim Vn ) 

\n— >oo / 

0(i.e., lim Vn has two different eigenvalues). 

n— >oo 

Then the system f[T^ has a basis of solutions itn^ and iti? with the 
following asymptotics as n ^ oo: 

= (^n[l+P'^/^1.2(^)]j (^2 + 0(1)), 
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where yUi and ^2 (Rg /ii < Re ^2) o^^e the eigenvalues of matrix V := 
lim Vn, "a^i and ~^ 2 ore the corresponding eigenvectors, fii{n) and 

/i2(^) are the eigenvalues of matrices Vn (chosen in a way such that 
Hiin) — )■ jii and fi2{n) — )■ /i2 as n ^ 00). 

Proof. The elhptic case of discrV^ < (a relatively straightforward one) 
is a special case of Janas-Moszynski Theorem proved in [9]. In the hy- 
perbolic case of discry > 0, the existence of the "smaller" solution 
^ is also guaranteed by the named theorem. So we need to prove 
only the existence of the second ("larger") solution itn^ (the solution 
corresponding to the eigenvalue fi2 with the largest real part) in hyper- 
bolic case. We emphasize that we do not give a new proof of the result 
from [9], but only add one extra (and rather simple) assertion to it. 

Let us reduce the situation to its simpler subcase, i.e., the system of 
the special form: 



(13) 



Un+l 
Vn+1 



For such a system, = V 



Vn 


-1 




+ Rr 



Ur. 

Vr, 



Remark 5. In fact we deal with a system of this type in Section 0, 
Step 3. 

The following statement holds: 
Lemma 1. Let 

• Pn be a positive for sufficiently large n sequence such that 

00 

J2Pn = +00. 

n=l 

• {Rn} ^ l^ be a 2 X 2 matrix sequence. 

Then there exists a solution to the system f[T^ of the form 



Un 
Vn 



?2 + 0(l), ^2 







as n 00. 



Proof. Without loss of generality one can assume that < < 1 for 
every n. Then for any two natural numbers ni < n2, 



'^2 

n 

n=ni 



I - i g ] + Rn 



n2 
n=ni 



I - 



Pn 





+ ll-Rnl 



"2 

< l[[l+\\Rj]<oo. 

n=ni 
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Therefore, every solution of the system ( |T3l) is bounded and there exists 
a universal constant C such that for any natural numbers < n2 and 
any solution m„ 



(14) 



Using the variation of parameters method one can rewrite the system 
(fT3|) in the following way: 





< c 




V ^"2 J 







(15) 



fc=no 

1 



u 



no 




n (i-pi) 

=k+l 

1 



Uk 
Vk 



The equivalence of the two systems (fT2!) and (fT5l) follows from elemen- 
tary calculations. Let us take 



u. 



then 



Note that H (1 

l=k+l 



no 



■'no 



n 

E 

k=no 



n a-pi) 

l=k+l 

1 



Vk 



Pi) — > as n — )■ oo for every /c due to properties 



of the sequence {p„}. Then since 



Un 



{Rn} G l^, one has by Weierstrass Theorem: 

Rk 



< C for n > uq and 




n i^-pi) 

l=k+l 

1 



fc=no 




1 



Rk 



Uk 
Vk 



Hence, 



Un 
Vr,. 



E 

fc=?10 




1 



Rk 



Uk 
Vk 



as n — )■ oo. The second component of the limit vector in the last 
expression is non-zero provided that tiq is sufficiently large (due to (fT^ 
and the second condition of Lemma [1]) . □ 
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Reduction of system ( |T2l) to the special case of system ( fT3l) can be 
done following the standard strategy [3], [H] by using the fact of D^- 
diagonalizability of every sequence of matrices with invertible limit 
[8] (also [9l Lemma 1.3], this is a discrete version of the result from 

my- 

Proposition 1. [8J Let {Vn} be a complex 2x2 matrix sequence such 
that {Vn} E and discr ( lim Vn) 7^ 0. Then the sequence {Vn] is 



D^-diagonalizable, i.e., there exists such a matrix sequence G 
with invertible limit that for sufficiently large values ofn, 



V = T 



j2i{n) 
fi2{n) 



T. 



-1 



Let us return to the proof of Theorem [21 By Proposition [T], for 
n sufficiently large the corresponding matrices Tn diagonalize Vn- To 
avoid tedious notations, without loss of generality we will do all the 
calculations starting with n = 1. An explicit calculation shows: 



T-Ul+PnVn + Rn)Tn 



where 
Further, 



Pn 



/ii(n) 

/i2(n) 
= I +Pn 



^ + T,^-\-lRnTn 

















/i2(n) 



Qni 



Tn+lRnTfi- 



Yl{I + PkVk + Rk 



k=i 



T 



n+1 



n 



I +Pk 



fll{k) 
/i2(fc) 



ll{l+Pkf^2m]Tn+llll 



where 



Pn ■= Pn- 



Kk=l >- 




/i2(n) - /ii(n) 

l+Pnf^2{n) 



1 Rn 



1 +Pn/^2(n) 



Qn- 



16 



SERGUEI NABOKO AND SERGEY SIMONOV 



The properties of the sequence {T„} guarantee that {Qn} £ as well 
as {Rn} G l^. Obviously pn > for large values of n, so the system 

satisfies all the conditions of Lemma [H Therefore it has a solution 
( ) = + 0(1). Let T := lim T„. One has: T~^2 = '^2 is an 

eigenvector of the matrix V corresponding to the eigenvalue 1x2- Then 
the system (JT2l) has a solution equal to 



\{[l+Pkl^2m]{'^2 + o{l))=: It^^ 



(2) 



which completes the proof. □ 
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